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I. ABSTRACT 

Intoraotlvo image data processing has greatly improved t^ho capability to extract 
scientific information from satellite imago data. Previous ■'>'ethoda of batch processing 
did not allow the scientist to interpret Intermediate results on a short leHoUgh time 
scale to effectively modify either the data or processes as re<)Uired for efficient data 
extraction. 

A primary Goddard Space Flight Center meteorological program objective is to develop 
methods to improve detection and prediction of severe storms using geostationary satellite 
mcasuroments. Interactive data processing systems can perform many tasks towards 
meeting this objective. 

The development of interactive image processing tcchniguos which have led to 
meteorological data processing systems at the Goddard Space Flight Center are discussed. 
Hardware and software functions of the existing Imoge Display and Manipulation System 
(IDAMS) and the special-purpose METPAK software package which has been developed for 
performing meteorology operations on IDAMS are described, A second-generation inter- 
active image processing system, the Atmospheric and Oceanographic Information Processing 
System (AOIPS) . which Is under development, is also outlined. These systems arc being 
designed to provide the necessary high throughput required for analysis of geostationary 
satellite data in the time domain. 


II. INTRODUCTION 

For over ten years, meteorologists at the Goddard space Flight Center (GSFO) have 
analyzed satellite data with computer assistance. The ability to modify and analyze 
data interactively is relatively new and offers to the user rapid step-by-step analysis 
with intermediate result presentation. Hardware and software systems discussed in this 
paper are specifically designed to interactively assist the meteorologist with the 
extraction and analysis of his data, 


III. GSFC METEOROLOGY PROGRAM OBJECTIVES 

Meteorological image data processing program objectives at Goddard Space Flight 
Center are centered on the anolysis of geostationary satellite data to obtain a quali- 
tative understanding of, and derive quantitative parameters for, the following: 1) the 

detection and prediction of severe storms, 2) the improvement of short-term and long- 
term global atmospheric numerical prediction models, and 3) the development of improved 
satellite systems. The driving requirement for a qualitative understanding is the 
multispectral display of a series of time-lapsed, enhanced Images of a rapidly developing 
severe local storm. The principal requirement of the quantitative parameter derivation 
is the measurement of cloud motions for the determination of wind velocity. Wind 
velocity measurements are essential for all types of analysis from the small-scale 
severe storm to the global atmospheric model. 


, 
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IV, DACKGHOUND 


In January o£ 1971 « an in-houBo dovalopud iiiiage procoBOlng oyatam was coroplotod 
ab Goddard Spaco Plight Contor. This aystom was doslgnatcd IDAMQ {Imago Display and 
Manipulation system) . IDAMS was dovolopcd to intoractivoly perform goomobrlc and 
radiomotric modification to digital imago data and bo display tho rosults on color 
and black-and-whlta tolovislon monitors, Opotabor Intoractlon with tho computer soft- 
ware and display hardware is achiovod through a computer graphics terminal, a 
joystick, and function control switches (boo Figure 1). Early in 1974, tho pobontial 
of using IDAMS for motoorological imago data processing was explored and tho decision 
mado to implement motoorological dato processing functions on IDAMS. Kathode used by 
universitios and other research organizations for processing of moteorological data 
wore evaluated prior to IDAMS software implementation. 

Two approoches to wind volocity (speed and direction) determination using tele- 
vision and computer systems hove been developed under NASA meteorology study contracts, 
one by Stanford Research Institute and the other by the University of Wisconsin, Doth 
methods are based on tho premise that wind vectors relate directly to change in 
cloud position over a period of time. Satellite image data rocoivod at known inter- 
vals provides a data base from which cloud sots are selected for tracking. 

The Stanford Research Institute system (SRI/MASA cloud study console) records 
photographic hard-copy imagery on an analog television "instant replay" disk recorder. 
Tima lapse viewing of the recorded imago scries, using conventional tolovislon monitors, 
enables the operator to select cloud sets which con be tracked to generate tho cloud 
motion wind osbimates. A cursor is positioned, by the operator, on each cloud to be 
tracked and the cursor position is transferred to punched paper tape for input to the 
computer. This process is repoatod on each successive image to establish cloud 
motion. Landmark data is assembled in a like manner and transferred to tho computer 
to be used as a geometric standard for imago data registration and geographic location. 
Cloud and landmark positions arc computer processed to produce an output tabulation 
and vector plot of the wind vectors, 

Tho University of Wisconsin system, designated HcIDAS (Man-Computer Interactive 
Data Access System) , is similar in function to the SRI wind vector generation; however, 
in this case, the processing computer is an integral part of tho system. This 
system approach offers considerably improved operator Interoctlon, system accuracy, 
and flexibility. Tho McIDAS input data base is from direct recorded satellite 
digital data. Digital data from approximately 30 SMS (Synchronous Meteorological 
Satellite) visible spectrum images (15,2SS 6-bit picture elements per line by 
14,568 lines on each image) and 30 Infrared images (3,822 8-bit picture elements per 
lino by 1,821 lines per image) are stored on a single 14-inch reel of magnetic tape. 
Image data are merged with orbit, landmark, and sensor related measurements in a pre- 
processing mode to establish tho geographic location of picture elements. 

Once geographic locations of the picture elements are known, precise positions of 
clouds selected to establish wind vectors are also known. Interactive cloud selection 
is followed by vector generation, which yields an image display on the television 
monitors of the original image data overlaid with arrows which represent the wind 
vectors. 

A cloud height preprocessing algorithm is being tested at the University of 
Wisconsin to determine the altitude of clouds used for wind Vector generation, thus 
providing wind vector data for several altitudes. The construction of three- 
dimensional wind vector sets will be possible using this approach. 

DEVELOPMENT OP FACILITIES 

A . IDAMS SYSTEM DEVELOP HPir't 

Design of the Image and Manipulation System (IDAMS) was initiated 

in early 1970 to provide an imaoe data manipulation capability primarily ror earth 
resources satellite data processing. Basic system guidelines required tho use of an 
existing computer manufactured by Control Data Corporation (CDC 3200) and the 
development of an interactive television terminal to be connected to the computer. 

The IDAMS interactive terminal was designed around a digital disk recording 
system which is used to refresh the television monitor screen 30 times per second. 



Appcoximatoly olx inllllon bita or throe tolovlaion images (500 x 700 x 5 bits) are 
stored In digital form on the dlaki Figure 2 is a simple block diagram of the VDAMS 
display terminal and computer ayatem. 

Image data is transferred from computer tope or disk through the comput-or Inter- 
face and o buffer memory to the IDAMS tolovlaion refresh disk. Transfer time is 13 
Boconda per imago. The throe stored images may bo displayed In true color, falso 
color, pseudo color, or black and white as the operator chooses. Digital inversion 
capabilitios allow o mixture of positive and negative image data for false color 
prosontation and image subtraction. A fade control is used to view two Images 
simultaneously. Image alternation for time lapse viewing and split soroon proeentatlons 
ore included in the terminal hardware. Selection of image alternation rates ond the 
sequence which the stored images are to be shown in time lapse is under either computer 
or manual control. A video disk recorder is available as the image storage media for 
time lapse viewing of up to 600 images. The video recorder is an analog system and 
therefore modification of image data must be performed prior to recording. An IDAHS 
hardware digital level slice function is used to divide the imago gray scale into three 
color display boundaries. The lower gray scale values are displayed in shadon of rod. 
the middle values in shades of green, and the upper levels in shades of blue. Two 
thumbwheel switches are used to select the dividing points between colors at any point 
on the 32-lovol gray &calc. 

Graphic overlay hardware functions include a crosshair or cursor and a variable 
sized rectangle. The position of the cursor and the rectangle's position, size and 
shape ate under direct computer control. The operator can position or modify rectangle 
size by either typing his input requirements or by. using the joystick and switches on 
the terminal console. The joystick and switches arc linked back to the computer so that 
their functions may be determined by software assignment. The IDAMS software package 
is very broad and encompasses most imago data manipulation operations being performed 
at image processing facilities around the country. Image display software is based 
on image reduction or zoom of input image data to television- sized images. Image data 
sots which are larger than approximately 500 X 700 picture elements are reduced for 
television monitor display by linear interpolation between adjacent elements. Smaller 
data sots are either displayed one-to-one or the zoom function is performed by repeating 
picture elements to .fill tho 500 X 700 dimension. Normal operations of the reduce ond 
zoom functions are as follows : 

1. Reduce large image data set for viewing 

2. Select an area of interest from the displayed image with tho graphic 
rectangle 

3. Zoom selected imago area to television screen dimensions and display 

In each case of reduce and zoom, original image data is used in the process to present 
the highest quality display and to preserve data integrity. 

Other IDAMS software functions may be used in a similar manner. An aroa Is 
operator selected from a reduced image for radiometric or geometric manipulation 
of the corresponding image data. The reverse of this approach is also a normal 
operation} that is, a television- sized image may be used for interactive manipulation 
and then, upon command, the same manipulation is applied to a larger imago data sot. 

Some of the significant general-purpose IDAMS software includes; linear and non- 
linear contrast modification; histogram development; image ratiolng; image subtraction: 
image data filtering; geometric rotation, projection, and correction. 

The general-purpose IDAMS radiometric and geometric image data manipulation 
software has been augmented with several special-purpose software systems. Those 
special software systems are discipline-related composites of existing imago 
manipulation software and special function software necessary to perform the user's 
operation. The IDAMS multispoctral classification and the meteorology cloud tracking 
wind vector software systems fall into this category, 

B. IDAMS METEOROLOGICAL DATA PROCESSING 

The special-purpose software package developed for performing meteorology operations 
on IDAMS is called MSTFAK. METPAK is an enhanced combination of modified McIDAS 
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cloud bracking and wind vector gcnorablon functions and IDAMS spccial-purpoBo Imago 
pcocosaing functions. The METPAK design is predicated on grouping of imago manipulation 
functions and bho oliroinotlon or minimization of housokceping raquirementst such as tapo 
assignmonts, file numbering, and other non-meteorology functions. 

The software system communicates with the operator by moans of menus and quostiona 
which require responses. Most rosponsoa ore in the form of pressing a single button to 
express operator decisions, 

METPAK functions for cloud trocking to gonorabo wind vectors are os follows > 

1, Establish latitude and longitude of selected clouds to accurately measure cloud 
position change from one imogo to another 

2, Compute wind vectors from the displacement of the selected clouds 

3, Display the wind vector field for evaluation and interactive modification 

4, Provide printer output for direct analysis and magnetic tapes of wind vectors 
for plotters and hard copy Imago generation systems 

The following paragraphs describe a typical cloud tracking - wind vector generation 
operation on IDAMS, 

Prior to a sosslon, the meteorologist selects a working imago data set by indicating 
data acquisition date, time, and the acquiring satellite. Data tapes are selected from 
the tape archive and sent to the IDAMS facility for processing. Image data sot size 
reduction to television display dimensions is either done as a non-prime time batch 
operation or may be accomplished using METPAK, 

The initial METPAK menu presented bo the operator, on the alphanumeric display, is 
shown in Figure 3, The Image data reduction sequence for television display is initiated 
when the aperotor typos bho number 1. The METPAK menu is replaced with bho image pro- 
cessing menu also shown in Figure 3, The image reduce function is selected by the oper- 
ator to proceed. The imago processing menu is replaced on the keyboard display with a 
series of instructions. The first Instruction is to have the computer operator mount 
the master data tape on a specific tape drive. The second instruction requests the 
operator bo name the reduced image. During each operation, informative statements are 
displayed to indicate to the user what is in process and where data is being stored. The 
reduced image data is automatically transferred to the IDAMS television refresh disk for 
display. 

The reduce function is repeated for three imago data sets which were acquired from 
the satcllibo at snlocted intervals. More data seta may be used for wind vector 
generation, however, the first image reductions are performed to allow the user to 
select landmarks and areas of meteorological interest. The throe images are viewed by 
the operator in a time lapse moving picture display to note position of viewable land, 
storms, and weather feature motion. The IDAMS graphic rectangle is used as a locator 
to outline the areas of interest for enlarge or zoom functions. 

The zoom function is first performed on areas the operator had selected to establish 
geographic earth reference. The operator selects the landmark register function from 
the METPAK menu to begin the landmark zoom sequence. The METPAK menu is replaced with 
the Landmark Registration menu from which the landmark zoom function is selected. 
Inquiries, such an zoom ratio, number of images to increase, and name of each image, are 
responded to by the operator and the processes are performed based on operator responses, 

Zoom ratios of either 1:1 or 4:1 are typical for locating landmarks. The 1:1 zoom 
ratio transfers one original data picture element to one television display clement. 

The 4:1 zoom transfers original data picture elements and lines four times for a magni- 
fied image display. 

After landmark zoom, the operator continues the landmark registration sequence by 
typing the number 3 (see Figure 3) . The operator positions the cursor over o known 
landmark on each zoomed image and assigns latitude and longitude coordinates to the 
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soloctod point or points, An accurate map facility must bo maintained for this purpose. 
A londmark dato file is occumulatod in tho computer as the operator locatoo landmarks 
on an imago sorios. When landmark soloction is complote, tho MtvTPAK master menu is 
again dlsployod. Tho operator types tho number 3 to initiate tho navigation process. 

The navigation program performs tho computations nocos<>ary to precisely locate 
earth coordinates of any clouds which will bo solocted £ai wind vector generation. 

Inputs to tho navigation programs ore satellite orbit pavamoters and landmark 
coordinates from tne landmark data file. The guality of, navigation Is tested and a 
navigation quolity report is prosontod on tho alphnnumorlc display. Poor navigation 
quality indicates improper landmark selection or incorrect orbit parameters. 

Corrections arc made by the operator, if roquirod, before proceeding to the cloud 
tracking function. 

Cloud tracking is tho foundation upon which wind vector generation from 
satellite data is based. Accurocy of wind vector sots generated by this method is 
based on the quality of ava^oble navigation data and the ability of the user to select 
clouds which move with tho wind. Tho operator typos tho number 4 to begin cloud 
tracking. 

The original reduced images are viewed in time lapse motion to aoloct an area of 
meteorological interest. Area selection and zoom are performed as previously 
described. Tho zoomed images are displayed in time lapse motion for selection of good 
cloud tracers which will bo used for wind vector generation. 

Two modes of cloud selection are possiblo using IDAMS METPAK. One mode is 
single-point tracking in which tho operator solects the center or an identifiable 
feature on a cloud with the cursor on three successive images. Point selection 
accuracy is enhanced by alternately switching between images to display the point 
selected on tho previous image and allowing tho operator to modify tho position of 
selected points on any of the images, viewing tho three images in time lapse motion 
with the selected points on each cloud displayed assures good single-point soloction. 
Tho second mode of cloud soloction is performed in a similar manner to single- point 
tracking, except thot a rectangle is used to outline solocted cloud areas on each of 
the three Images. i:loud displacements ore computed by best fit cross correlation 
between the three cloud image data sots. This method requires loss effort by the 
operator and is more precise than single-point tracking but requires mere computer 
time. 


As cloud selection progresses in either mode, the operator is kept aware of 
selected clouds by the cursor momentarily flashing at each point previously selected. 
This function is necessary to eliminate duplicate selection and to insure complete 
cloud selection coverage over the displayed area. Upon completion of cloud selection, 
wind vector programs are initiated by typing the number 5. 

Wind vectors are calculated by applying tho previously determined navigation 
coordinate transform to tho picture coordinates of a selected cloud sot and establishing 
true earth displacement (distance and direction) of the clouds. Cloud displacement 
is divided by the time interval between images to calculate cloud velocity. Arrows 
are generated to present wind vector information to the user with arrow length 
proportional to wind speed. Arrow position and pointing direction relate wind 
direction as established from a particular cloud set. 

Wind vectors ar« the output product of the IDAMS METPAK and arc presented in 
several forms for use by the meteorologist. The display mode is used to observe wind 
flow and to edit out any obvious erroneous vectors. Display options include: vector 

field display, or vector display overlaid on selected images. Time lapse motion 
viewing of the vector overlaid Image sets may bo used for dynamic analysis of atmos- 
pheric motion. 
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Another wind vector output product is a computer printout o£ the wind vector 
information. The printout lists dotO) tlno, latitudoi and longitude of each cloud 
Bolootod and the wind apeod and direction from each cloud sot. 

Magnetic tapes of the vector and image overlaid vector data are another output 
product. These tapes are used to produce hard-copy products . such asi ink plot 
vectors for map overlay and photographic outputs of the vector sots overlaid on tho 
original imago series. Tho choice of output products is made by tho operator from a 
menu similar to other HETPAK operations. 

A utility program set has boon included as part of tho METPAK sy8t6m to allow 
an operator to manipulate imago data prior to. or during a processing session. The 
utility programs are intended to remove some of tho rigidity imposed by tho user* 
oriented METPak programs, use of tho utility programs is monitored to determine now 
system requirements in order to implement now METPAK and utility operations. 

Data managoment programs are used throughout tho HETPAK software system to track 
imago operations and stored results. Other routines within the data management 
programs allow an operator to stop and restart from any point in tho METPAK operation. 
This provision is made to ollmlnato duplication of effort from session to session. 

Tho operator may terminate an operation at tho end a aossion and reinitiate ot the 
same point in his next nossion. saved and restart d.'tta is preserved on digital 
magnetic tape for re-entry. 

IDAMS is being used as a technique development system for meteorology and earth 
resources applications. The METPAK software package has enabled the mctcorolOuV 
users to easily extract scientific data interactively with minimal system knowledge, 

Tho system must be limited to technique dovolopmont, however, since an operational 
system requires much higher data throughput than is possible on IDAMS. IDAMS will 
bo used for continued meteorology data processing technique development through 
early 1976. An operational system, the Atmospheric and Oceanographic Information 
Processing System (AOIPS) , is currently under development, AOIPS will perform 
motoocolcgy and earth resources data processing as well as new technique development 
on a much larger scale thon the current IDAMS effort. Tho ostimated starting date of 
operation is November 1975. 

C. TUB ATMOSPHERIC AND OCEANOGRAPlIlO INFOKMATION PR0CES5INS SVSTEH (.AOIPS) 

AOIPS is a dual terminal multiprocessor image data pcccosslng system with 
much greater data manipulation and throughput capocity thon tho IDAMS facility. 
Meteorology and earth resources applications users share AOIPS for their interoctivo 
data manipulation analysis functions. 

Central processing and control ore performed by two Digital Equipment Corporation 
PDP 11 computers interconnected as multlprcccssors. Either computer may use any 
of the peripheral equipments in both foreground and background modes of operation. 

System managoment software is used to assign equipment priorities for optimum data 
throughput. The system block diagram shown in Figure 4 illustrates interconnection 
of tho AOIPS terminals, computers, and related equipment. One of the terminals is 
scheduled for delivery in August of 1975. Tho other terminal is to bo delivered in 
early 1976. Several other items shown in the block diagram are also under development 
with a‘..ticipatcd delivery in late 1975 or early 1976. The present configuration hos 
a modified General Electric Imago-100 serving as the interactive terminal. 

The Imago-100 system has been motiifiod to allow ..equonce viewing of an image 
BoricB and to incorporate computer control of all hardware functions necessary for 
image manipulation. D!>'.gltal imago data la stored in tho Imago-100 terminal for tele- 
vision screen refreshing and hardware manipulation. Five solid-state memories provide 
storage for five television sized images. Each memory can store one eight-bit image 
with dimensions of 512 picture elements per lino and 512 lines per imago. The Imoge-lOO 
terminal is directly connected to one of the POP 11 computers. Data linkage between 
computers is accomplished by interconnecting the computer data buses and certain control 
functions. 

Data tape archives ore maintained to provide users with an accosalble data 
base. Standard 9-track digital magnetic tapes and hlgh-dcnslty digital tapes are 
stored in close proximity to the AOIPS facility and In the Goddard central tape 


Btorago. Oporations poceoimol aro roaponslblo for acquiring data tapes prior to a user 
soBsion and roturning thorn to tho appropriate storage after use. 

The AOIPS hlgh'dansity digital tape recorder will bo used for imago data input. A 
largo data comprosBion ratio is achieved by using hlgh-donaity tapes rather than stand- 
ard 9-track computer tapes. Image data from 250 9-track tapes may be stored on one hlgh- 
donsity tape. Savings of tape, storage space, and data acquisition time are consider- 
able. 

Tolcvision-sizod images generated during an operating session are stored for rapid 
recall on the two largo computer disk systems. Accumulation of over COO digital tclo- 
vision-slzod images is possible , Data management software allocates storage area to 
each user for later recall of his imago data. Larger imago data sets may bo stored 
should a User want to manipulate more than tho viewed Imago area. 

AOIPS system software has boon under development since early 1975. Tho motoorology 
software package equivalent to MEXPAK will bo ready by December 1975. Data management 
software is also under development to provide data acquisition and storage operations on 
a large scale. 

Upon delivery of the second Interactive terminal, tho system will bo fully linked 
for multitcrminal and multiprocessor operation. The facility will provide the users with 
continuous 24-hour service. Preparatory operations will be performed at night to allow 
the Users full utilization of their allocated time for application-oriented data 
manipulations. A staff of system specialists will optimize data flow based on user needs 
to further facilitate system use for scientific analysis, 

VI. SUMMARY AND FUTURE WORK 

Interactive image data manipulation has allowed effective analysis of meteorological 
satellite data in the time domain which would otherwise be impossible, Future efforts 
will concentrate on tho following i 

1. Refinement of data processing techniques and meteorological analysis 
functions 

2. The development of new meteorological analysis functions 

3. The addition of a large on-line computer to allow tho retrieval of data 
from large data banks and to perform analysis functions on a larger scale 

4. The development of a coal-time facility for severe storm forecasting 
research 
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Figure 2. Image Display and Manipulation System (IDAMS) 
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Figure 3. METPAK Master Menu and Sub-Menu 
Examples 
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Figure 4. AOIPS Block Diagram 





















